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THE OPTICAL LUMINOSITY FUNCTION OF VIPJALIZED SYSTEMS 
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ABSTRACT 

We determine the optical luminosity function of virialized systems over the full range of density en- 
hancements, from single galaxies to clusters of galaxies. The analysis is based on galaxy systems identified 
from the Nearby Optical Galaxy (NOG) sample, which is the largest, all-sky catalog of objectively- 
identified bound objects presently available. 

We find that the B-band luminosity function of systems is insensitive to the choice of the group-finding 
algorithms and is well described, over the absolute-magnitude range —24.5 < M + 51og/i75 < —18.5, by a 
Schechter function with a s = -1.4±0.03, M* - 5 log /175 = -23.1 ±0.06 and (j>* = 4.8 x 10" 4 h% 5 Mpc" 3 
or by a double power law: <j> p \(L s ) oc £-i-45±o.07 f Qr ^ < ^ an( j 0(£ s ) oc £-2.35±o.i5 f or ^ > ^ 

with L p i = 8.5 x 10 10 /i 75 2 L Q , corresponding to M s — 5\oghjs = —21.85. The characteristic luminosity of 
virialized systems, L p \, is ~ 3 times brighter than that (L* al ) of the luminosity function of NOG galaxies. 
Our results show that half of the luminosity of the universe is generated in systems with L s < 2.9L* al 
and that 10% of the overall luminosity density is supplied by systems with L s > 30L* a! . 

We find a significant environmental dependence in the luminosity function of systems, in the sense 
that overdense regions, as measured on scales of 5ft, -1 Mpc, preferentially host brighter, and presumably 
more massive, virialized systems. 

Subject headings: cosmology: large-scale structure of the universe — galaxies: clusters: general — 
galaxies: halos — galaxies: luminosity function 



1. INTRODUCTION 

From a theoretical perspective, dark matter halos are 
the basic building blocks of the structure in the Universe. 
Observationally, halos must be identified with virialized 
galaxy "systems" , ranging from single dwarf galaxies to 
rich clusters of galaxies. The missing link between the two 
is the mass-to-light ratio, which in general is a function 
of mass or luminosity. As a first step toward understand- 
ing this relationship, in this paper we determine the lu- 
minosity function, not of individual galaxies, but of entire 
virialized systems. This can be compared with the lumi- 
nosity functions theoretically predicted by semi-analytic 
models of galaxy formation and with other halo statistics 
(Marinoni & Hudson 2001), such as the mass function. 

While the luminosity function (LF) of cluster and group 
galaxies has been extensively studied in literature (Iovino 
et al. 1993; Marzke & da Costa 1997; Lumsden et al. 1997; 
Gaidos 1997; Rauzy, Adami, & Mazure 1998; Bromley et 
al. 1998; Ramella et al. 1999; Marinoni et al. 1999), only 
few detailed studies have examined the LF of galaxy sys- 
tems. Gott & Turner (1977, GT) used data for 103 groups 
identified in two dimensions (i.e. from the angular coordi- 
nates of galaxies only) to determine the functional shape of 



the observed group LF. Bahcall (1979), combining the GT 
data with the Abell (1958) clusters, attempted to extrap- 
olate the observed luminosity distribution into the rich- 
cluster domain. These early studies, however, suffered 
from a number of weaknesses: groups identified only in 
projection on the sky, small-number statistics, and uncer- 
tain completeness and corrections. These weaknesses were 
largely overcome by MFW, who determined the LF of sys- 
tems using 163 groups (with > 3 members) identified in 
the CfAl redshift survey (Davis et al. 1982; Huchra et al. 
1983), the best data set available at that time. The qual- 
ity and quantity of the data have improved considerably 
since those studies, allowing a more robust determination 
of the LF of virialized systems (hereafter VSLF). 

Our VSLF is based on the Nearby Optical Galaxy 
(NOG) catalog (Marinoni 2001), a nearly all-sky, com- 
plete, magnitude-limited sample of ~ 7000 bright and 
nearby galaxies. In order to assess the effects of group- 
ing algorithms, we will analyze three catalogs of groups, 
all based on NOG but extracted using different algo- 
rithms and selection criteria (Giuricin et al. 2000). When 
grouped, NOG contains ^2800 isolated galaxies and ~ 
1100 systems (~600 binaries and ~500 groups with at least 
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three members). This group catalog combines a represen- 
tative volume (cz < 6000 km s -1 , nearly all-sky) of the 
Universe with a high comoving density of galaxies to al- 
low excellent statistics in the group regime. 

This paper is the fifth in a series (Marinoni et al. 1998, 
Paper I; Marinoni et al. 1999, Paper II; Giuricin ct al. 
2000, Paper III; Giuricin ct al. 2001, Paper IV;) in which 
we investigate the properties of the large-scale structures 
as traced by the NOG sample. The outline of our paper 
is as follows: in §2 we review and summarize the iden- 
tification procedures of NOG galaxy systems and discuss 
some specific properties of NOG groups. In §3 we describe 
our method for calculating the VSLF and apply it to our 
sample. Results are summarized in §4. 

Throughout this paper, the Hubble constant is taken to 
be 75 /175 km s _1 Mpc -1 and the recession velocities czlg 
are evaluated in the Local Group rest frame. 

2. THE NOG GROUPS: IDENTIFICATION AND SPATIAL 
DISTRIBUTION 

The NOG sample is a statistically controlled, distance- 
limited (czlg < 6000km s _1 ) and magnitude- limited (B < 
14) complete sample of more than 7000 optical galaxies. 
The sample covers 2/3 (8.27 sr) of the sky (\b\ > 20°), 
a volume of 1.41 x 10 6 h^ 5 Mpc -3 and has a redshift com- 
pleteness of 98% (see Marinoni 2001). In this paper we use 
the small extension of the NOG (for a total of 7232 galax- 
ies) containing 156 additional galaxies which have rough 
estimates of magnitudes (see Paper III for details). Note 
that, in contrast to most previous studies of large-scale 
structure, the magnitudes used in this paper are homog- 
enized total blue magnitudes, given in the standard sys- 
tem of the RC3 catalog (de Vaucouleurs et al. 1991) and 
fully corrected for Galactic and internal extinction and K- 
dimming. 

Redshift maps present two principal distortions: i) 
small-scale perturbations caused by random velocities in 
clusters of galaxies which produce the so-called "Fingers of 
God" i.e. a radial stretching in galaxy maps pointed at the 
observer and ii) large-scale perturbations caused by large 
overdensities leading to coherent bulk motions. For the 
purposes of this paper, the latter effect of large-scale pe- 
culiar velocities are not severe; these can be easily treated 
via the flow models derived in Paper I. 

We corrected for the Finger-of-God effect in Paper III, 
where galaxy systems were identified by means of different 
objective, group-finding algorithms i.e. the widely-used hi- 
erarchical (H) and percolation (P) (or friends- of -friends) 
algorithms. The P algorithm (Huchra & Geller 1982) iden- 
tifies as members of galaxy aggregations the galaxies which 
have their transverse separation Dij < Do ■ R and line-of- 
sight velocity difference czij < cz ■ R 1 where R is a scal- 
ing parameter which takes into account the decrease of 
the magnitude range of the luminosity function sampled 
at increasing distance. In Paper III, two variants of the 
percolation method were used to identify groups. In one 
variant (hereafter denoted as "scaled" friends-of- friends 
or P2, following Paper III), typical values for the linking 
parameters, at the median depth of the sample (~ 4000 
km s _1 ), are 434 km s _1 in velocity and 0.89 /17s -1 Mpc in 
transverse separation, and both link parameters were then 
suitably scaled with distance to give a density threshold 



of — = 80. In the other variant (hereafter denoted sim- 
ply as percolation or PI), only the transverse separation 
link parameter was scaled with distance, while the veloc- 
ity link parameter was fixed at the value of 350 km s _1 . 
Given the limited range of redshift encompassed by the 
NOG, this choice was used to approximate a slow scaling 
of the velocity link parameter with distance, as suggested 
by cosmological N-body simulations (e.g., Nolthenius & 
White 1987). 

The hierarchical clustering method introduced by 
Materne (1978) and revised by Tully (1980) is an algo- 
rithm which generates a hierarchical sequence of systems 
organized by some affinity parameter. In paper III, the au- 
thors adopted pl = 8 x 10 9 L Q Mpc~ 3 (corresponding to a 
luminosity density contrast = 45) as the limiting lu- 
minosity density parameter used for cutting the hierarchy 
and defining groups. 

For both algorithms, the adopted values of the group- 
selection parameters were chosen after a search in the pa- 
rameter space guided by numerical simulations in order 
to obtain realistic and homogeneous catalogs of groups. 
Since all three catalogs are based on the same galaxy sam- 
ple, this allows us to investigate systematic effects related 
to the choice of grouping algorithm. 

Most of the NOG galaxies (~60%) are found to be mem- 
bers of galaxy binaries (which comprise ~15% of galaxies) 
or groups with at least three members (~45% of galaxies). 
About 40% of the galaxies are left ungrouped (isolated 
galaxies) . 

In Paper III, the similarity between the catalogs of 
groups was demonstrated using the fraction of members 
in common to groups of different catalogs. In this paper, 
we also show that the three catalogs of groups show sim- 
ilar statistical distributions and trace similar large-scale 
structures. Figure 1 shows the distributions of P2 groups 
on the sky. The smoothed group number density contrast 
at a given point r (in units of velocity) has been computed 
by summing the contribution from all the i groups of our 
catalog: 
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are the mean number density, smoothing and selection 
functions respectively, and where mn m = 12 is the limiting 
apparent magnitude of completeness of the group catalog 
(see §3.2). Smoothing is performed in three dimensions 
with a variable Gaussian smoothing length R s given by 
the mean inter-group separation ranging from 300 to 670 
km s _1 . 




Fig. 1. — NOG galaxy systems reconstructed using the scaled 
friends of friends (P2) algorithm are shown in equal-area Aitoff 
projection on the sky using Galactic coordinates. The region de- 
void of galaxies is the zone of avoidance (\b\ < 20°). The thick line 
is drawn at the celestial equator, <5 = 0°. Galaxy systems are indi- 
cated by filled circles, where the symbol size scales with the number 
of members of each groups. Also shown is the smoothed distribution 
of systems on shells at rcdshift 1000 km s~ 1 (uppcr), 3000 km s _1 
(center), 5000 km s _1 (lower). Smoothing is performed in 3D with a 
variable Gaussian smoothing length given by the mean group sepa- 
ration (which is given by 309, 480 and 622 km s -1 respectively). The 
heavy contour denotes the mean density. Dashed contours represent 
densities below the mean (spaced with A(^) = 0.2 intervals). Re- 
gions in excess of the mean are in gray scafe with contours spaced 
with A(^) = 0.4 intervals. 




Fig. 2. — NOG isolated galaxies, i.e. galaxies left ungroupcd after 
running the scaled friends of friends algorithm (P2), are shown in 
equal-area Aitoff projection on the sky using Galactic coordinates. 
The region devoid of galaxies is the zone of avoidance (|6| < 20°). 
The thick line is drawn at the celestial equator, <5 = 0°. Also 
shown is the smoothed distribution of systems on shells at redshift 
1000 km s~ 1 (upper), 3000 km s _1 (center), 5000 km s — 1 (lower). 
Smoothing is performed in 3D with a variable Gaussian smoothing 
length given by the mean group separation (which is given by 309, 
480 and 622 km s _1 respectively). The heavy contour denotes the 
mean density. Dashed and solid contours represent densities below 
and above the mean and are spaced with A(^) = 1/3 intervals. 
Regions in excess of the mean are in gray scale. 
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Fig. 3. — Cone diagrams of the distribution of NOG galaxies in 
the Perseus-Pisces region: the distance coordinate is the recession 
velocity (in the Local Group frame), angular coordinates is Galac- 
tic longitude (1). Upper:) distribution of galaxies before grouping; 
Lower Left:) distribution of objects after the application of the hi- 
erarchical (H) group-finding method; Lower Right:) distribution of 
objects after the application of the friends of friends (PI) group- 
finding algorithm. The same dot scaling as in Figure 1 is used to 
represent groups with different number of members. 




Fig. 4. — As in Figure 2 but for the Hydra-Centaurus region. 
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A comparison between Fig. 1 and Fig. 2 shows that 
groups delineate the same large-scale structures as do 
galaxies but with a different density contrast, i.e. groups 
are "biased" with respect to galaxies. In a separate paper 
(paper IV), we presented a detailed analysis of the clus- 
tering of NOG groups. 

In Figures 3 and 4, we use cone diagrams in order to 
illustrate graphically the performance of two of our recon- 
struction algorithms (the H and PI methods) in two inter- 
esting regions of the sky, the Perseus-Pisces and Hydra- 
Centaurus superclusters. These plots show that the dif- 
ferent catalogs of groups have consistent spatial distribu- 
tion. The top cone shows the ungrouped distribution of 
galaxies in redshift-space. Virialized systems show up as 
"Fingers of God" which are elongated along the observer's 
line of sight. Most groups as well as isolated galaxies are 
recovered by both methods. Moreover, these plots demon- 
strate how well our optimized grouping algorithms perform 
across a wide range of density enhancements and across the 
different large-scale structures in which the galaxies are 
embedded (i.e. orthogonal and parallel to the observer's 
line of sight). 
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Fig. 5. — Upper: distribution of 1093 hierarchical binaries and 
groups as a function of distance in the northern and southern 
(shaded histograms) Galactic hemisphere. Lower: fraction of galax- 
ies in hierarchical groups as a function of distance across the overall 
NOG volume. Data are binned in 500 km s _1 distance intervals. 

Figure 5 illustrates that the fraction of NOG galaxies 
that are grouped is roughly a constant, independent of 
distance, although, due to large-scale clustering, the num- 
ber of systems is an increasing function of distance in the 
southern galactic hemisphere and a decreasing one in the 
north. This indicates that, although we lose an increas- 
ing number of galaxies and groups with distance, due to 
the magnitude-limited sample, our group reconstruction 
methods do not suffer from any serious distance-dependent 
bias. Fig. 5 shows the results for H groups, but the other 
catalogs of groups give similar results. 



3. THE LUMINOSITY FUNCTION OF VIRIALIZED SYSTEMS 

The NOG groups are drawn from a magnitude-limited 
galaxy sample. In order to calculate the VSLF, we regard 
all systems as aggregates of galaxies much in the same way 
as galaxies are aggregates of stars. However, the crucial 
difference is that while galaxy fluxes can be directly mea- 
sured with some device down to a limiting surface bright- 
ness, the flux of a system has to be inferred from mem- 
bers above a given apparent magnitude. When calculat- 
ing the LF of galaxies, we must first correct for missing 
stars, which fall below the surface brightness limit, and 
then correct for the missing galaxies which fall below the 
magnitude limit of the catalog. In the same way, when 
constructing the LF of systems, we must first correct for 
missing galaxies (which fall below the magnitude limit of 
the survey) and then correct for missing systems. 

We assume here that every system contains galaxies 
drawn from a universal galaxy LF, and thus has a large 
number of faint members. Consequently, even a single 
galaxy should be considered as the sole member observed 
above the apparent magnitude limit of a group, and thus 
requires a correction for unseen luminosity. By applying 
a luminosity correction to all systems, we obtain a "viri- 
alized systems" catalog. We stress that the result is fun- 
damentally different from the group catalog from which 
it is derived. The group catalog suffers from the well- 
recognized problem that two groups with N visible ele- 
ments identified (from a magnitude-limited galaxy sample) 
at two different distances correspond to potential wells of 
different masses. The systems catalog is more physical in 
the sense that it corresponds to halos with total absolute 
luminosity L. 

To overcome these difficulties and estimate the lumi- 
nosity function of systems, we must 1) assume a model 
describing how the missed galaxies are clustered around 
"seeds" defined by NOG systems, 2) check the photomet- 
ric completeness of the resulting catalog in order to en- 
sure that its luminosity distribution is representative of 
the group population and not only of our specific sam- 
ple; 3) use a statistical estimator which is independent of 
density fluctuations because we know that groups tend to 
cluster at least as strongly as do galaxies (see Fig. 1). 



3.1. The Luminosity Selection Function 

We begin by introducing the statistical method used for 
estimation of the total luminosity L s of systems. 

The NOG galaxy LF, corrected for Malmquist bias and 
calculated after correcting for clustering is described by 
a Schechter-type function (Schcchtcr 1976) with the fol- 
lowing shape parameters: a ga \ = —1.10 ± 0.03, M* al — 
51og/i 75 = -20.61 ±0.04 (Marinoni 2001). 

Marinoni (2001) also showed that the LF of galaxies 
in groups is marginally consistent with the field galaxy 
LF. Therefore, we assume that the observed luminosity in 
each bound system is a random realization of that por- 
tion of the universal LF <f>(L) which lies above the local 
magnitude limit at the group distance, evaluated from 
the median redshift of its galaxy members. Specifically, 
Lum = 10°- 4 ( M o- M >- )L Q , where M Q = +5.48 mag and 
Mim = mn m - 5 log(cz is ) -15 + 5 log h 75 . 
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Fig. 6. — Upper: Spaenhauer diagram (Spacnhauer 1978) show- 
ing magnitude versus redshift for groups with N > 2 members. 
The absolute magnitude of the systems is the total one, i.e. cor- 
rected for unseen members. The solid line is the expected pre- 
diction for groups without any correction applied and corresponds 
to the following apparent variation in absolute magnitude M a 



51og{c2 LG > - 15 + 51og/i 75 (m s ' 



13.25). The dot- 



ted line is the visibility function of the NOG galaxy catalog and 
corresponds to an apparent magnitude cutoff = 14. Lower: 

The additive luminosity corrections per galaxy, in units of L*y as 
given by our luminosity selection function is compared with the 
MFW method, using the magnitude system adopted in the NOG 
and the magnitude cut-off Tig™ = 14. The hatched area shows how 
our luminosity selection function varies when it is computed using 
a power-law approximation at the faint end (M fainter than —16) 
with an exponent in the range from —1.1 to —1.8. 

In order to recover the total number of galaxies in a sys- 
tem, the observed number density of galaxies is weighted 
by wn{t) = 1/S(r), where S(r) is the selection function of 
a magnitude-limited sample (cf. eq. 4), yielding a corrected 
number density n(r) = n oos (r)wN(r). For luminosity, we 
can write the total expected luminosity of a system as 



L s (r) = L ohs (r) + L(r) = L oha (r)w L (r) 



(5) 



where L denotes the unseen luminosity from galaxies below 
the magnitude limit , and wl{t) is the luminosity- density 
weighting function (see also Gourgoulhon, Chamaraux, & 
Fouque 1992) 



/ °° L(j>(L)dL 



L4>(L)dL 



(6) 



As r goes to zero, the luminosity weighting function 
iiJi(r) reduces to 1 as desired. With this weighting scheme, 
we obtain an estimate of the absolute luminosity of groups 
which is independent of the magnitude limit of the catalog 
in which systems arc identified, is independent of density 
fluctuations (i.e. independent of the LF normalization fac- 
tor (f)*) and is less sensitive to the faint-end uncertainties 
and environmental dependencies of the galaxy LF. 



There are other ways in which the total group light 
can be recovered. For example, MFW assumed the ob- 
served number of galaxies in each group to be a ran- 
dom sample of the observed portion of the LF and 

L = 



additive correction term of the form 
^L^dL/d^^dL). 



used an 
Nobs Jo u 

Our correction and that of MFW are compared in Fig- 
ure 6 where we plot the fractional expected correction in 
L* al units as a function of distance for a single galaxy with 
luminosity L* al . Note that we do not expect these two cor- 
rections to be identical, since our correction is applied to 
the observed luminosity, whereas that of MFW is applied 
to the observed number. 

We have tested our correction by moving nearby groups 
to progressively greater distances. However, at large dis- 
tances our approach has the advantage that it depends 
very weakly on the adopted Schechter parameters because, 
at the faint end, the correction does not weight using the 
number density of dwarf galaxies, but using their luminos- 
ity density. (Note that our correction does not diverge 
using as the lower limit of integration.) Thus, even 
if the galaxy function suffers some environmental depen- 
dence, i.e. becoming steeper at the faint end in clusters 
(Marzke & da Costa 1997), this correction is insensitive 
to this change. This can be seen explicitly in Figure 6, 
where we show the effect on our correction when we mod- 
ify the faint-end behavior of the Schechter galaxy luminos- 
ity function by adding a power-law term with an exponent 
equal to —1.8 for luminosities fainter than M = —16 (see 
Zucca et al. 1997). The correction to the total system lu- 
minosity is small (< 20%) and is not a strong function of 
distance. As noted above, we find a shallow slope (—1.1) 
for NOG galaxies and find no evidence for a large differ- 
ence in the NOG B-band LF parameters for field galaxies 
versus galaxies in groups and clusters. Thus we can take 
the above correction as a conservative upper limit on such 
systematic effects. 

3.2. Sample Completeness 

We have applied eq. 6 to all systems, including isolated 
galaxies. This yields three "virialized system" (VS) cata- 
logs of objects which, depending on the catalog of groups 
used, contain ~ 4000 objects (~ 2800 systems with one 
observed member and ~1100 groups with at least 2 mem- 
bers). For example, the hierarchical VS sample is com- 
posed by 2828 isolated galaxies and 1093 groups (with 4404 
members) for a total of 3921 "systems". There are 691 bi- 
nary systems, 297 systems with 3 < N < 5 members and 
195 systems with N > 5 members. 

The large number of VS considered is not by itself a 
guarantee that our catalogs are photometrically complete. 
Instead, it is necessary to analyze the scaling of the number 
of systems as a function of the reconstructed magnitude 
and consider as suitable systems for the LF estimate only 
those for which the observed logarithmic scaling follows 
roughly the expected Euclidean N s cx I0 6ms behavior 
(no cosmological corrections are needed within the NOG 
volume). Each VS subsample (singles, binaries, etc.) has 
been analyzed with this method. As an example, in Figure 
7, we plot integral and differential (independent bins) log- 
arithmic counts for the sample of 195 hierarchical systems 
containing N > 5 members. 
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Fig. 7. — Integral (left) and differential (right) counts (binned 
in 0.4 magnitude intervals) for NOG groups with iV ga i > 5 mem- 
bers. Apparent magnitudes have been corrected with the luminosity 
selection function. Bars represent ±1(7 Poisson errors. Solid lines 
indicate the expected prediction in a volume with an Euclidean ge- 
ometry and a homogeneous distribution of systems. 

The Euclidean growth rate is followed over a range of 3 
magnitudes and breaks down at B ~ 11.1 mag. Roughly 
the same behavior (but with different limits of complete- 
ness) is observed in different subsamplcs of the same VS 
catalog or different VS catalogs. The magnitude limits 
of completeness and the number of objects which meet 
these criteria are reported in Table 1 for the case of the 
H groups. Although this selection reduces the number of 
systems available for the analysis, systematic errors due to 
incompleteness are avoided. 



3.3. Results 

In recovering the LF shape and normalization for vari- 
ous group subsamplcs, we have used the estimation tech- 
nique adopted in Paper II. In particular, the shape has 
been derived through Turner's (1979) method and the nor- 
malization has been calculated using the relation <j>* = 

n s / J M ^(j)(M s )dM s with M c = mf m - 5 log 500 - 15 + 

51og/i75, and where rn£ m is the limiting apparent mag- 
nitude of completeness of each subsample and n s has 
been determined using the minimum variance estimator 
of Davis & Huchra (1982). 

Our estimate of the VSLF is shown in Figure 8 for the H 
and PI VS samples. We have found that a Schechter func- 
tion with parameters a s = —1.4 ± 0.03, M* — 51ogft, 7 5 = 
-23.1 ± 0.06 and </>* = 4.8 x 10~ 4 h% 5 Mpc~ 3 provides a 
good fit to data over a broad range of absolute luminos- 
ity (-24.5 < M s - 5logh 75 < -18.5). In Figure 8 we 
also show the joint distribution of errors in a s and M* as 
derived from the x 2 matrices of the least squares fit. 
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Fig. 8. — The virialized systems luminosity function (differential 
form) as derived using NOG hierarchical (H) (solid circles together 
with ±1(7 error bars) and friends of friends (PI) (open triangles) 
group catalogs. The best fitting Schechter LF function (solid line) 
and the 1, 2, 3 a confidence ellipses are also shown. The solid tri- 
angle represents in the (a B — M*) space the best-fitting Schechter 
parameters of the friends of friends LF. 

Figure 8 also shows that different catalogs of groups (H 
and PI) do not give significantly different LFs. Moreover, 
since the P2 groups (obtained by varying the velocity link- 
ing parameter) give a LF which is very similar to the LF 
obtained from the PI groups in which the velocity link 
parameter is kept fixed, we argue that while other group 
properties (e.g. velocity dispersion) considerably depend 
on the values chosen for this parameter (e.g. Trasarti- 
Battistoni 1998), the LF is rather stable and only weakly 
sensitive to different choices. Besides we can use redshifts 
as distance indicators and obtain a VSLF which is not dif- 
ferent from a LF that would be obtained in a pseudo-real 
space analysis. For the sake of simplicity, in the following 
we shall present results obtained using the H groups. 

While the LF is rather insensitive to the group-finding 
algorithm, it is somewhat sensitive to the choice of density 
threshold. If thresholds are set much higher than the ones 
used in this paper, this yields a larger proportion of field 
galaxies and, hence, to a steeper faint-end of the LF, in 
agreement with the results of MFW. Marinoni & Hudson 
(2001) show that our adopted density thresholds lead to 
systems which are virialized. 

In Figure 9 we show that the VSLF can be economically 
described in terms of a smoothly decreasing double-power 
law of the following form 



(f>(L s )dL s 



^ te) 



-1.45±0.07 
-2.35±0.15 



d(^) if L s <L pl 
'd(±-) if L s >L pl 



(7) 



10 




-24 -22 -20 -18 

M s -5 log h 75 



Fig. 9. — The VS luminosity function is plotted (solid circles) and 
compared to the CfA VSLF of MFW (squares) and to the NOG 
galaxy LF (dotted line). Solid lines represent a double power law fit 
to data. The 1, 2, 3 a confidence ellipses of the best fitting Schechter 
LF are also shown. The solid triangle represents in the (ce s — M*) 
space the best-fitting Schechter parameters of the CfA group LF as 
given by MFW (after having transformed Zwicky magnitudes into 
the B system used in the NOG). 

where L p i = 8.5 x 10 10 h^L ( 7 ) , corresponding to M p \ = 
-21.85 and </> pl = 6.5 x 10~ 4 h% 5 Mpc~ 3 . In contrast with 
the galaxy LF, which has a sharp break at L* the VSLF 
shows a more slowly-varying behavior. In the limit of low 
luminosities, the VSLF approaches the NOG galaxy LF 
as expected, since these "systems" are single field dwarf 
galaxies. 

The global mean density of the VS distribution is 
n <~ 2.2 x 10~ 2 /i 3 5 Mpc~ 3 for systems with a magni- 
tude brighter than M s — 5\ogh 75 = —15 and n <~ 6.2 x 
10~ 3 /i 3 5 Mpc -3 for systems brighter than M s — 5 log /175 = 
-18. 

In Figure 9, we also show the good agreement (< 2a) be- 
tween the CfA groups LF and the NOG VSLF. In order to 
compare our results (which use the corrected B total mag- 
nitude in the RC3 system) to the results of MFW (which 
are based on the Zwicky magnitude system, B z ), we have 
made the following transformations: B — B z = —0.35 mag 
(Auman et al. 1989) and -E>( corr ) — B = —0.25 mag (for the 
internal absorption, see Paper II). 

The two volumes are largely independent, with NOG 
being shallower and wider in area than the CfAl volume 
studied by MFW. The agreement between the two VSLFs 
is reassuring because the CfA groups have been identified 
using a variant of the P algorithm with different selection 
parameters 




M s 5 log h 75 

Fig. 10. — Composite group LF obtained by summing the con- 
tribution from halos hosting visible systems of different richness. 
Points represent the VSLF. The corresponding 1 a error ellipses are 
plotted in the inset. 

and scaling laws, as well as a different luminosity cor- 
rection and a different method for estimating the LF. 
Moreover, if we consider the GT VSLF, we see that their 

bright-end slope has a power-law behavior (L s ) with 
essentially the same slope of our LF. However, due the 
lack of small systems with low surface density enhance- 
ment with respect to the background, their data show a 
faint end which is flatter than our LF and that of MFW. 
Our results are in marginal agreement with the LF de- 
rived by Bahcall (1979) who finds a significantly flatter 
slope (~ 3a) at the faint end. 

In Figure 10, we have estimated the LF for different 
range of group richness from isolated galaxies to groups 
with N> 5 members. The best-fit parameters are given in 
Table 1. While field galaxies dominate the low-luminosity 
portion of the VSLF, richer systems contribute in a pro- 
gressive way to the bright end of the VS LF. We also show 
the total LF defined as the sum over the LFs pertaining 
to different samples. 

The remarkable flatness of the luminosity-weighted 
VSLF was first noted by GT. A significant contribution 
to the total luminosity density of the universe is made 
by both small and large systems (see Figure 11). The 
characteristic luminosity of the NOG galaxy LF is L* al = 

2.7 x 10 10 /i 75 2 L Q We find that after correcting for the un- 
seen luminosity, systems with total luminosity L s < L* al 
make up 25% of the total luminosity density pl of the uni- 
verse, systems with L s < 5L* al and L s < 10£* a] contribute 
by 58% and 75% respectively, while a substantial fraction 
of luminosity density (10%) is contributed by large systems 
with more than 30L* ,. 
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Fig. 11. — Fraction of the total luminosity density from systems 
with L < L s as derived from the VSLF. L*^ is the characteristic 
luminosity of the NOG galaxy LF (M*, - 5 log /i 75 = -20.61) while 
L B (M b ) is the absolute luminosity of galaxy systems. 



3.4. Environmental Effects on the Luminosity Function 
of Virialized Systems 

The mass function of dark matter halos is expected to 
show a dependence on the large-scale environment, as pre- 
dicted by e.g. Mo & White (1996) using an extension of the 
Press-Schecter formalism, and confirmed in N-body cold 
dark matter simulations (see also Lemson & Kauffmann 
1999). The sense of the predicted trend is that dark mat- 
ter halos in overdense regions should be biased toward the 
high masses than those forming in lower density regions. 
Note that the extended Press-Schecter theory makes pre- 
dictions for the clustering of virialized dark matter halos, 
which are typically groups or clusters rather than individ- 
ual galaxies. 

Many workers have looked for such an effect by studying 
the luminosities of galaxies as a function of environment, 
with conflicting results (e.g., Valotto et al. 1997; Lopez- 
Cruz et al. 1997; Bromley et al. 1998; Ramella et al. 1999; 
Zabludoff & Mulchaey 2000; Christlein 2000; Balogh et al. 
2001). Part of this may be due to different bandpasses used 
in different studies, the environmental dependence being 
strong in the i-C-band (e.g. Balogh et al 2001) and weak 
in the B-band. Marinoni (2001) compared NOG galaxies 
in regions of high density contrast to those in low density 
contrast, and also compared isolated galaxies to members 
of systems. He found that the £>-band luminosity func- 
tion of galaxies as well as the luminosity functions of early 
and late morphological types are consistent with no depen- 
dence on environment, although a weak trend (differences 
in M* of a few tenths of a magnitude) is also permitted by 
the data. 

It is clear that at the very bright end, more luminous 
galaxies are more strongly clustered (Norberg et al. 2001), 
and hence one should expect some dependence of the LF 
on environment. The trend found by Norberg et al, while 
highly significant, is rather weak: the biasing parameter 
6/6* = 0.85 + 0.15-L/.L*, where 6* is the biasing parameter 
of an L* galaxy. 




-24 -22 -20 -18 

M s 5 log h 75 



Fig. 12. — The VS luminosity function for hierarchical systems in 
regions of different density. The environmental density parameter 
^ for each halo has been weighted using a Gaussian filter (eq. 3) 
with a fixed smoothing length R s = 530 km s _1 which corresponds 
to the NOG mean inter-group separation. The inset shows the 1, 
2, 3 cr contours for the joint distribution of errors of the Schechter 
parameters (a B ,M*) of systems in low density regions. The solid 
triangle represents in the (a B — M*) space the best-fitting Schechter 
parameters of the LF of systems in high density regions. 

Thus at a magnitude brighter than L*, the relative bias- 
ing is only 1.22. This leads to a rather weak effect of when 
comparing, for example, Schecter function parameters in 
overdense and underdense regions. 

In contrast to the situation for individual galaxies, we 
expect that the total luminosities of virialized systems are 
closely related to the dark matter content of the virial- 
ized halo. Systems should therefore be expected to dis- 
play the environmental dependence predicted by the ex- 
tended Press-Schecter theory. We can test this predic- 
tion as follows. We assign to each system an environmen- 
tal density parameter (^-) computed by smoothing the 
data with a Gaussian filter (eq. 3) of fixed smoothing 
length R s = 530 km s _1 (which corresponds to the mean 
NOG inter-group separation). Systems which lie in under- 
dense regions (— < 0; 361 objects) and overdense regions 
(— > 0; 946 objects) are differently distributed in lumi- 
nosity. Note that the system in question is excluded from 
the smoothed density estimate. 

Figure 12 shows that the hierarchical VSLF is sensitive 
to the large-scale environmental density: the luminosity 
distribution of halos which inhabit low density regions is 
biased toward faint magnitudes. 

Both the distributions are well described in term of a 
Schechter function with parameters [-1.8,-23.3] and [-1.2, 
-22.7] respectively in low and high density regions. The 
difference in Shcchter function parameters is significant at 
the 3<7 level. A signal of the same type is also present, al- 
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beit at a marginal (1.5(7 ) significance level when we smooth 
the halo density field on a larger scale (R s = 1000km s _1 ). 

Thus, in contrast to the situation for galaxies, the VSLF 
is clearly dependent on the large-scale environment. This 
effect is not the result of the rather weak dependence of 
B-band galaxy LFs on environment. More luminous, and 
presumably more massive, systems are more abundant in 
high-density regions and hence are more strongly clus- 
tered, in agreement with the group correlation function re- 
sults of Paper IV. One astrophysical implication of this re- 
sult is that the higher relative biasing of early-type galax- 
ies is then a byproduct of the environmental dependence 
of virialized halos coupled with the morphology-density 
relation. 

4. SUMMARY AND CONCLUSIONS 

The traditional approach has been to classify non-linear 
structures as being either galaxies, groups or clusters. In 
contrast to either galaxies or rich clusters of galaxies, how- 
ever, groups have received little attention in the literature, 
despite the fact that they these virialized objects contain 
most of the luminosity, and hence presumably most of the 
mass, in the universe. Our approach has been to treat all 
virialized "systems", from isolated single galaxies to rich 
clusters of galaxies, as a continuum, and to consider the 
total optical luminosity of these systems. 

We have analyzed the NOG group catalog (Paper III) in 
order to correct the total B-band luminosities of these sys- 
tems for the luminosity of galaxies below the magnitude 
limit of the sample and for the completeness of subsam- 
ples of the group catalog, and have used this to measure 
the luminosity function of virialized systems. We verify 
that the luminosity function is insensitive the choice of 
the group-finding algorithm. It is also insensitive to the 



value of the velocity link parameter adopted in the perco- 
lation algorithm, whereas it is somewhat sensitive (at the 
faint end) to the density contrast at which groups are de- 
fined. Our luminosity function is in good agreement with 
that of MFW, who used a quite different approach to its 
determination and relied on groups identified in the CfAl 
survey. 

The luminosity function of virialized systems is well de- 
scribed by a double-power law, (p(L s )dL s cx L~ 1A5±0 - 07 dL 
for L s <L* and (j){L a )dL s oc L- 2 - 35±0A5 dL s for L s > L* 
with L* = 8.5 x 10 10 /if 5 2 L Q , over a broad range of abso- 
lute luminosity (-24.5 < M s - 5 log h 75 < -18.5). Our 
results indicate that 25%, 50%, and 75% of the luminosity 
of the universe is generated in systems with L s < -£<* al , 
L s < 2.9L* al and L s < 10L* gal , respectively; 10% of the 
overall luminosity density is supplied by systems with 
L s > 30L* al . 

Finally, we find that the luminosity function of systems 
depends on the large-scale environment, in the sense that 
halos distributed in low-density regions host preferentially 
low-luminosity systems, while there is an excess of bright 
virialized systems in high-density regions. 
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Table 1 

Luminosity function parameters for various samples of systems 





All systems 


Isolated 
galaxy 
"systems" 


Binary 
systems 


Groups 3 < N g <5 


Groups N g >5 




12. 


12. 


12. 


11.5 


11.1 


N s 


1307 


377 


450 


247 


187 


v/v max 


0.47 


0.61 


0.54 


0.49 


0.31 


a 


-1.40 ±0.03 


-1.7 ±0.12 


-1.2 ±0.12 


-1.2 ±0.15 


-0.8 ±0.20 


M* B -5logh 75 


-23.1 ±0.06 


-21.7 ±0.30 


-22.2 ±0.22 


-22.6 ±0.27 


-23.0 ±0.28 


0*(1CT 4 ^ 5 Mpc~ 3 ) 


4.8 


3.9 


2.9 


1.6 


2.7 


p L (l0 8 L Q h 3 75 Mpc" 3 ) 


2.0 


0.86 


0.39 


0.32 


0.59 



